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Production of new superheavy Z=108-114 nuclei with 238 U, 244 Pu and 248,250 Cm targets 

Zhao-Qing FengQ Gen-Ming Jin, and Jun-Qing Li 

Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, People's Republic of China 

(Dated: December 21, 2009) 

Within the framework of the dinuclear system (DNS) model, production cross sections of new 
superheavy nuclei with charged numbers Z=108-114 are analyzed systematically. Possible combina- 
tions based on the actinide nuclides 238 U, 244 Pu and 248 ' 250 Cm with the optimal excitation energies 
and evaporation channels are pointed out to synthesize new isotopes which lie between the nuclides 
produced in the cold fusion and the 48 Ca induced fusion reactions experimentally, which are feasible 
to be constructed experimentally. It is found that the production cross sections of superheavy nuclei 
decrease drastically with the charged numbers of compound nuclei. Larger mass asymmetries of the 
entrance channels enhance the cross sections in 2n-5n channels. 

PACS number(s): 25.70.Jj, 24.10.-i, 25.60.Pj 



The synthesis of superheavy nuclei (SHN) is motivated 
with respect to search the "island of stability" which is 
predicted theoretically, and has obtained much experi- 
mental research with fusion-evaporation reactions [3, 0] • 
Neutron-deficient SHN with charged numbers Z=107- 
112 were synthesized in cold fusion reactions with the 
208 Pb and 209 Bi targets for the first time and investi- 
gated at GSI (Darmstadt, Germany) with the heavy- ion 
accelerator UNILAC and the SHIP separator [l|, Q . Re- 
cently, experiments on the synthesis of element 113 in the 
70 Zn+ 209 Bi reaction have been performed successfully 
at RIKEN (Tokyo, Japan) Q|. More neutron-rich SHN 
with Z=113-116, 118 were assigned at FLNR in Dubna 
(Russia) with the double magic nucleus 48 Ca bombarding 
the actinide nuclei [2j, [5|, l6|. New heavy isotopes 259 Db 
and 265 Bh were also synthesized at HIRFL in Lanzhou 
(China) 0]. New SHN between the isotopes of the cold 
fusion and the 48 Ca induced reactions are of importance 
not only for investigating the structure of SHN such as 
influence of shell effect on stability of SHN etc, and also 
as a stepstone for further synthesizing and identifying 
heavier superheavy nuclei. 

Here we use a dinuclear system (DNS) model [1, 
in which the nucleon transfer is coupled to the relative 
motion by solving a set of microscopically derived master 
equations by distinguishing protons and neutrons, and 
a barrier distribution in the capture and fusion process 
of two colliding nuclei is introduced in the model. In 
the DNS model, the evaporation residue cross section 
is expressed as a sum over partial waves with angular 
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Here, T(E c , m ,, J) is the transmission probability of the 
two colliding nuclei overcoming the Coulomb barrier in 
the entrance channel to form the DNS. The Pcn is the 
probability that the system will evolve from a touching 
configuration to the formation of compound nucleus in 
competition with the quasi-fission of the DNS and the 
fission of heavy fragment. The last term is the survival 
probability of the formed compound nucleus, which can 
be estimated with the statistical evaporation model by 
considering the competition between neutron evapora- 
tion and fission [8j. 

Within the concept of the DNS, the fusion probabil- 
ity was also calculated by using the multidimensional 
Kramers-type expression to get the fusion and quasifis- 
sion rate by Adamian et al. [ill ]. In order to describe 
the fusion dynamics as a diffusion process along proton 
and neutron degrees of freedom, the fusion probability 
is obtained by solving a set of master equations numer- 
ically in the potential energy surface of the DNS. The 
time evolution of the distribution probability function 
P(Z\, Ni, Ei, t) for fragment 1 with proton number Z\ 
and neutron number N\ and with excitation energy E\ 
is described by the following master equations, 



dP{Zl, j t 1,El,t) =]T W Zl , Ni;ZltNl (t) [d Zl >Nl P(Z[,N 1 ,E[,t)-d zi!Nl P(Z 1 ,N 1 ,E 1 ,t)]+J2w Zl ,N 1 -,z 1 ,N l (*) 

[dz^P^N^E^t) - d ZltN ^P(Zx, Ni, Ex, t)] - [A,/(e(t)) + A /is (0(i))] P{Z\,N\,E\,t){2) 
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Here W Zlt N 1 - z ^N 1 {W Zi ,Ni;Z x ,n'^) is the mean transition probability from the channel (Zi,Ni,Ei) to (Z^,Ni, E[) 
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(or (Z 1 ,N 1 ,E 1 ) to (Z ll N[,E' 1 )), and d Zl . Nl denotes the 
microscopic dimension corresponding to the macroscopic 
state (Z X) N X ,E X ). The sum is taken over all possible 
proton and neutron numbers that fragment Z[ , N[ may 
take, but only one nucleon transfer is considered in the 
model with the relation Z' x = Z x ± 1 and N{ = N x ± 1. 
The excitation energy E x is determined by the dissipation 
energy from the relative motion and the potential energy 
surface of the DNS. The motion of nucleons in the inter- 
acting potential is governed by the single-particle Hamil- 
tonian @, Q . The evolution of the DNS along the vari- 
able R leads to the quasi- fission of the DNS. The quasi- 
fission rate A q f and the fission rate Aj is of heavy frag- 
ment are estimated with the one-dimensional Kramers 
formula [l,[lo|. 

In the relaxation process of the relative motion, the 
DNS will be excited by the dissipation of the relative 
kinetic energy. The local excitation energy is determined 
by the excitation energy of the composite system and 
the potential energy surface of the DNS. The potential 
energy surface (PES) of the DNS is given by 

U({a}) =B(Z X ,N X )+B(Z 2 ,N 2 ) - [B{Z, N) + V r c f (J)] 
+V({a}) (3) 

with Z x + Z 2 = Z and N x + N 2 = N. Here 
the symbol {a} denotes the sign of the quantities 
Zi,N x , Z 2l N 2 ; J, R; 0i,0 2 , 0i,9 2 . The B(Z ll N l )(i = 
1,2) and B(Z,N) are the negative binding energies of 
the fragment (Zi,Ni) and the compound nucleus (Z,N), 
respectively, which are calculated from the liquid drop 
model, in which the shell and the pairing corrections are 
included reasonably. The Vf?J* is the rotation energy of 
the compound nucleus. The (3i represent the quadrupole 
deformations of the two fragments. The 9{ denote the 
angles between the collision orientations and the sym- 
metry axes of deformed nuclei. The interaction poten- 
tial between fragment (Z X> N X ) and (Z 2 , N 2 ) includes the 
nuclear, Coulomb and centrifugal parts, the details are 
given in Ref. [9(. In the calculation, the distance R be- 
tween the centers of the two fragments is chosen to be the 
value which gives the minimum of the interaction poten- 
tial, in which the DNS is considered to be formed. So 
the PES depends on the proton and neutron numbers of 
the fragments. 

The formation probability of the compound nucleus at 
the Coulomb barrier B and for the angular momentum 
J is given by [1, Q 

Zbg Nbg 

PcN(E c . m .,J,B)= J2 Yl P(Zi,N x ,E x , nnt ). (4) 

Z 1= l Ni=l 

The interaction time Tint in the dissipation process of two 
colliding partners is dependent on the incident energy 
i? c m . and the quantities J and B. We obtain the fusion 



probability as 

PcN(E c . m „, J) — J f(B)P CN {E c . m .,J,B)dB, (5) 

where the barrier distribution function is taken as an 
asymmetric Gaussian form. 

Neutron-deficient SHN with charged numbers Z=107- 
113 were synthesized successfully in the cold fusion re- 
actions. The evaporation residues was observed by the 
consecutive a decay until to take place the spontaneous 
fission of the known nuclides, in which the fusion dynam- 
ics and the structure properties of the compound nucleus 
have a strong influence in the production of SHN. Re- 
cently more neutron-rich and heavier SHN with charged 
numbers Z=113-116, 118 were produced in the fusion- 
evaporation reactions of 48 Ca bombarding actinide tar- 
gets. Superheavy residues were also identified by the 
consecutive a decay, unfortunately to spontaneous fis- 
sion of unknown nuclides. Neutron-rich projectile-target 
combinations are necessary to be chosen so that super- 
heavy residues approach the " island of stability" with the 
doubly magic shell closure beyond 208 Pb at the position 
of protons Z=114-126 and neutrons N=184. New SHN 
between the isotopes of the cold fusion and the 48 Ca in- 
duced reactions are of importance for the structure stud- 
ies themselves and also as daughter nuclides for identify- 
ing heavier SHN in the future. 

The excitation energy of compound nucleus is obtained 
by Eq N = E c rn , + Q, where £ c , m , is the incident energy 
in the center-of-mass system. The Q value is given by 
Q — AMp + AMt — AMc, and the corresponding mass 
defects are taken from Ref. [T3| for projectile, target and 
compound nucleus, respectively. Usually, neutron- rich 
projectiles are used to synthesize SHN experimentally, 
such as 64 Ni and 70 Zn in the cold fusion reactions, which 
can enhance the survival probability W sur in Eq.(l) of 
the excited compound nucleus owing to smaller neutron 
separation energy. Within the framework of the DNS 
model, we calculated the evaporation residue cross sec- 
tions of superheavy element Mt based on the actinide 
targets 248 ' 250 Cm, 244 Pu and 238 U with the neutron-rich 
projectiles 27 A1, 31 P and 37 C1 as shown in Fig. 1. One 
can see that the 3n channel in the reactions 27 Al+ 248 Cm 
and 37 C1+ 238 U, and the 4n and 5n channels in the sys- 
tem 27 Al+ 250 Cm have the larger cross sections in the 
production of SHN 272 Mt and 273 Mt. Superheavy ele- 
ment Ds(Z=110) was successfully synthesized in the cold 
fusion reactions [H, UH • The production of the SHN de- 
pends on the isotopic combinations of projectiles and 
targets. Calculations were performed for the reactions 

30 Si+ 248,250 Cm) 36 S+ 244p u and 40 Ar+ 238TJ tQ produce 

superheavy element Ds as shown in Fig. 2. Combina- 
tion with 248 Cm has the larger cross section in the 4n 
channel than the isotope 250 Cm due to the larger value 
of survival probability. The 4n channels in the systems 

30 Si+ 248,250 Cm &nd 40 Ar+ 238TJ &nd the 3n channe l m 
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FIG. 1. Evaporation residue excitation functions in the production of isotopes of superheavy element Mt in the reactions 
27 Al+ 248 ' 250 Cm, 31 P+ 244 Pu and 37 C1+ 238 U. 
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FIG. 2. Calculated production cross sections for the reactions 30 Si+ 248 ' 250 Cm, 36 S+ 244 Pu and 40 Ar+ 23S U to produce superheavy 
element Ds. 



the reaction 30 Si+ 248 Cm are feasible in the synthesis of 
new SHN 274_276 Ds. These combinations can be chosen 
in experimental preparation with the present facilities. 

In the DNS model, the isotopic trends are mainly de- 
termined by both the fusion and survival probabilities. 
When the neutron number of the projectile is increasing, 
the DNS gets more symmetrical and the fusion probabil- 
ity decreases if the DNS does not consist of more stable 
nuclei due to a higher inner fusion barrier. A smaller neu- 
tron separation energy and a larger shell correction lead 
to a larger survival probability. The compound nucleus 
with closed neutron shells has larger shell correction en- 
ergy and neutron separation energy. The cross section 
decreases rapidly in the production of the isotopes of 
Rg(Z=lll). Optical channels are the 4n evaporation for 



the systems 31 P+ 248 > 250 Cm, 37 Cl+ 244 Pu and 41 K+ 238 U 
to produce 275 > 277 Rg a s shown in Fig. 3. Superheavy 
element Z=112 is more difficult to be produced in the 
selected systems. Shown in Fig. 5 gives that the possible 
way is the 4n channel in the reaction 36 S+ 250 Cm, but the 
cross section is still smaller than the system 48 Ca+ 238 U 
fioj ] although the larger mass asymmetry. 

The productions of superheavy element Z=113 were 
successfully performed in the cold fusion reaction 

70 Zn+ 209 Bi [| and alsQ in the ho< . fuf , ion 48 Ca+ 237 Np 

[14j | with the cross section less than 1 pb. We calculated 
the evaporation residue excitation functions for the reac- 
tions 37 Cl+ 248 ' 250 Cm, 41 K+ 244 Pu and 45 Sc+ 238 U. The 
results show that the 3n and 4n channels in the sys- 
tems 37 Cl-f 248,250 Cm have larger cross sections and are 
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FIG. 3. The same as in Fig. 1, but for the reactions 31 P+ 248 ' 250 Cm, 37 Cl+ 244 Pu and 41 K+ 23S U to produce superheavy element 
Rg. 




FIG. 4. Comparison of the calculated production cross sections in the reactions 36 S+ 248 ' 250 Cm and 40 Ar+ 244 Pu to synthesize 
superheavy element Z=112. 



possible to synthesize new isotopes 281 - 284 H3 i n exper- 
imentally. Superheavy element Z=114 is difficulty to be 
produced from our calculations for the selected systems 
because of the smaller cross sections with less than 0.1 
pb for all systems. We list the maximal production cross 
sections and the corresponding excitation energies in the 
brackets calculated by using the DNS model in Table 1. 
These selected systems and evaporation channels are fea- 
sible to produce new isotopes between the cold fusion and 
the 48 Ca induced fusion reactions. 

In summary, we systematically investigated the pro- 
duction of superheavy residues in the fusion-evaporation 
reactions using the DNS model, in which the nucleon 
transfer leading to the formation of superheavy com- 
pound nucleus is described by a set of microscopically 



derived master equations distinguishing the proton and 
neutron transfer that are coupled to the dissipation of 
relative motion energy and angular momentum. The pro- 
duction of new isotopes between the gap of the cold fu- 
sion and the 48 Ca induced fusion reactions are discussed 
for selected systems. Optimal evaporation channels and 
excitation energies corresponding to the maximal cross 
sections are stated and discussed systematically. 

We would like to thank Prof. Werner Scheid for care- 
fully reading the manuscript. This work was supported 
by the National Natural Science Foundation of China 
under Grant Nos. 10805061 and 10775061, the spe- 
cial foundation of the president fund, the west doctoral 
project of Chinese Academy of Sciences, and major state 
basic research development program under Grant No. 
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TABLE I. Comparisons of calculated maximal evaporation residue cross sections and optimal excitation energies (in bracket) 
in 2n-5n channels. 



Reactions 



at R (pb) (E* CN ) 



t (pb) (E* CN ) 



2b Mg+ 248 Cm 
26 Mg+ 250 Cm 
30 Si+ 244 Pu 
36 S_)_ 238 U 

27 A1+ 248 Cm 
27 A1+ 250 Cm 
31p + 244p u 
37 C1+ 238tj 

30 Si+ 248 Cm 
30 Si+ 250 Cm 
36 S+ 244 Pu 

40 Ar+ 238 u 

31 P+ 248 Cm 
31 P+ 250 Cm 

37 a+ 244p u 

41 K+ 238 U 
36 S+ 248 Cm 
36 S+ 250 Cm 

40 Ar+ 244p u 
37 a+ 248 Cm 
37 a+ 250 Cm 
41 K+ 244p u 

45 Sc+ 238 U 

40 Ar+ 248 Cm 

40 Ar+ 250 Cm 
48 Ti+ 238 U 
50 Ti+ 238 U 



2.50 (40 MeV) 
1.11 (41 MeV) 
0.46 (42 MeV) 
0.21 (37 MeV) 
0.31 (44 MeV) 
0.12 (46 MeV) 

4.71 xlO" 2 (47 MeV) 
0.16 (38 MeV) 
0.34 (39 MeV) 
0.11 (41 MeV) 

3.87xl0" 2 (36 MeV) 

0.26 (32 MeV) 
4.11 xlO" 2 (43 MeV) 
9.91 xlO -3 (47 MeV) 
4.01 xlO" 2 (36 MeV) 
1.96xl0" 2 (32 MeV) 
4.31 xlO" 2 (33 MeV) 
1.76xl0" 2 (36 MeV) 
5.79 xlO" 3 (31 MeV) 
5.81 xlO" 2 (31 MeV) 
2.08 xlO" 2 (35 MeV) 
1.11 xlO" 2 (29 MeV) 

1.72 xlO" 2 (27 MeV) 
6.98 xlO" 3 (26 MeV) 
2.77xl0" 3 (29 MeV) 
4.51 xlO" 2 (24 MeV) 
5.11 xlO" 2 (23 MeV) 



26.2 (40 MeV) 
10.57 (41 MeV) 
5.09 (43 MeV) 
1.96 (38 MeV) 
27.83 (44 MeV) 
1.64 (46 MeV) 
4.25 (47 MeV) 
13.31 (38 MeV) 
1.72 (39 MeV) 
0.42 (42 MeV) 
0.101 (38 MeV) 
0.55 (36 MeV) 
0.31 (44 MeV) 
5.49 xlO" 2 (48 MeV) 

0.11 (38 MeV) 
8.67xl0" 2 (35 MeV) 
7.64 xlO" 2 (36 MeV) 
8.25 xlO" 2 (37 MeV) 
9.48xl0" 3 (36 MeV) 
0.26 (35 MeV) 
0.21 (36 MeV) 
4.22 xlO" 2 (34 MeV) 
1.99xl0" 2 (35 MeV) 
2.21 xlO" 2 (33 MeV) 
2.11 xlO" 2 (33 MeV) 
1.37xl0" 2 (32 MeV) 
2.18xl0" 2 (31 MeV) 



719.1 (42 MeV) 

185.2 (42 MeV) 
185.1 (44 MeV) 
42.97 (42 MeV) 
3.59 (47 MeV) 

24.31 (46 MeV) 
0.87 (50 MeV) 
0.67 (44 MeV) 

65.32 (43 MeV) 
3.54 (44 MeV) 
0.61 (41 MeV) 
2.10 (42 MeV) 
1.85 (45 MeV) 
0.41 (51 MeV) 
0.33 (42 MeV) 
0.21 (41 MeV) 

7.02xl0 -2 (43 MeV) 

0.24 (42 MeV) 
9.84xl0 -3 (44 MeV) 
0.195 (42 MeV) 
0.594 (41 MeV) 
3.31 xlO" 2 (42 MeV) 
2.32xl0 -3 (45 MeV) 
3.5xl0" 2 (41 MeV) 
7.96 xlO" 2 (40 MeV) 
5.81 xlO" 3 (42 MeV) 
1.07xl0" 2 (40 MeV) 



1.23 (51 MeV) 
108.5 (45 MeV) 
0.72 (51 MeV) 
0.11 (52 MeV) 
1.34 (51 MeV) 
97.44 (49 MeV) 
0.52 (53 MeV) 
0.17 (50 MeV) 
1.22xl0" 2 (56 MeV) 
0.93 (48 MeV) 
0.12 (48 MeV) 
0.45 (53 MeV) 
0.69 (50 MeV) 
0.25 (52 MeV) 
7.15xl0 -2 (49 MeV) 
3.77xl0" 2 (49 MeV) 
2.55xl0 -3 (54 MeV) 
1.47xl0" 2 (51 MeV) 
5.02xl0 -4 (56 MeV) 
1.05xl0" 2 (52 MeV) 
7.06 xlO" 2 (49 MeV) 
2.3xl0" 3 (53 MeV) 
1.92xl0" 4 (57 MeV) 
2.12xl0" 3 (51 MeV) 
9.69xl0~ 3 (48 MeV) 
1.71 xlO" 4 (54 MeV) 
4.11xl0" 4 (50 MeV) 



2007CB815000. 



* fengzhq@impcas.ac.cn 
[1] S. Hofmann and G. Miinzenberg, Rev. Mod. Phys. 72, 

733 (2000); S. Hofmann, Rep. Prog. Phys. 61, 639 (1998). 
[2] Yu. Ts. Oganessian, J. Phys. G 34, R165 (2007). 
[3] G. Miinzenberg, J. Phys. G 25, 717 (1999). 
[4] K. Morita, K. Morimoto, D. Kaji et al, J. Phys. Soc. 

Jpn. 73, 2593 (2004). 
[5] Yu. Ts. Oganessian, V. K. Utyonkov, Yu. V. Lobanov et 

al, Phys. Rev. C 69, 021601(R) (2004). 
[6] Yu. Ts. Oganessian, V. K. Utyonkov, Yu. V. Lobanov et 

al, Phys. Rev. C 74, 044602 (2006). 
[7] Z. G. Gan, Z. Qin, H. M. Fan et al, Eur. Phys. J. A10, 



21 (2001); Z. G. Gan, J. S. Guo, X. L. Wu et al, Eur. 

Phys. J. A20, 385 (2004). 
[8] Z. Q. Feng, G. M. Jin, F. Fu and J. Q. Li, Nucl. Phys. 

A771, 50 (2006). 
[9] Z. Q. Feng, G. M. Jin, J. Q. Li and W. Scheid, Phys. 

Rev. C 76, 044606 (2007). 
[10] Z. Q. Feng, G. M. Jin, J. Q. Li and W. Scheid, Nucl. 

Phys. A816, 33 (2009). 
[11] G. G. Adamian, N. V. Antonenko, W. Scheid et al, Nucl. 

Phys. A627, 361 (1997); Nucl. Phys. A633, 409 (1998). 
[12] P. Moller et al, At. Data Nucl. Data Tables 59, 185 

(1995). 

[13] S. Hofmann, V. Ninov, F.P. Hefiberger et al, Z. Phys. A 
350, 277 (1995); Z. Phys. A 350, 281 (1995). 

[14] Yu. Ts. Oganessian, V. K. Utyonkov, Yu. V. Lobanov et 
al, Phys. Rev. C 76, 011601(R) (2007). 



